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Abstract. - The behavior of weakly deformed drops on non wetting surfaces is usually described 
using linear models. We show that these simple pictures cannot account for measurements of the 
dynamics of droplets that oscillate or bounce on super-hydrophobic substrates. We demonstrate 
that several peculiar experimental observations observed in previous works can be understood 
through a logarithmic correction of the linear model. 
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l Water impacts on surfaces are one the most common 
physical phenomena observed in every day life, with a 
wide range of applications, going from agriculture to inkjet 
printers. Recent progresses in experimental and theoret- 
ical methods have yielded a better understanding of this 
problem PQH]. In particular, the development of novel 
superhydrophobic surfaces has allowed the study of these 
phenomena in well controlled situations [2J0] , where some 
of the most challenging effects, and most notably the line 
contact dynamics and hysteresis [5], can safely be ne- 
glected. 

In this Letter, we present a description of the dynam- 
ics of a sessile droplet in the regime of small deforma- 
tion. Contrasting with the commonly accepted scenario, 
we show that the droplet behaves as a non-linear spring 
whose stiffness decreases logarithmically with the defor- 
mation. Our model is based on an analogy with the clas- 
sical Hertz contact for an elastic sphere and leads to ac- 
curate quantitative comparisons with existing data on su- 
perhydrophobic surfaces. We suggest that the analysis 
can be extended to the oscillations of droplets sitting on 
non- wetting substrates. 

Consider a droplet of radius R, surface tension a and 
density p at the surface of a non-wetting substrate. For 
R <C K , where = yojpg is the capillary length 
and g is gravity, then the sag e the drop (see Fig. [2]) 
is given by a balance between a linear capillary restoring 



force ~ as and the weight ~ pgR yielding the scaling 
relation e ~ k 2 R 3 0(7]. 

We measured the dimcnsionless sag ej^R? for water 
droplets coated with lycopodium grains (a way to generate 
super hydrophobicity [5]), and the results are presented in 
figure [1] As seen on this semi-log plot, this ratio is not 
observed to be independent on the drop radius for Rn < 1, 
but it weakly (logarithmically) diverges, suggesting rather 
a scaling 



s~ ^R A |ln(«i?)| 



(1) 



This singular scaling was first pointed out in a study 
of the static properties of wet foams [9], and it can be 
justified by calculating the shape of the drop in the regime 
of weak deformation using an expansion very similar to 
that of [TD1[TT]. Laplace equation relating the pressure 
discontinuity to the mean curvature C of the liquid/air 
interface reads in spherical coordinates (Fig. O 



pgr cos 6 = Pq + 2aC, 



(2) 



where Pq is the atmospherical pressure and Pq is an in- 
tegration constant for the hydrostatic equations. In the 
limit of small deformations, we write r = R + (, with 
C <C R, and we obtain |12j for an axisymmetric drop 
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(2 + A e ) C = AP + pgR cos ( 



(3) 
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Fig. 1: Color online. Measurement of the sag of the top of a 
liquid marble obtained by deposition of lycopodium spores at 
the surface of a water droplet. Black dashed line : theoretical 
calculation using the numerical integration of the exact Laplace 
law [2] Red solid line : asymptotic expansion (Eq. [7]) in the 
weakly deformed regime. This asymptotic expansion is valid 
only in the regime kR < 1 where the drop is weakly deformed. 



where APq = Pq — Pq — 2<t/R and Ag is the angular Lapla- 
cian given by 



sin(0) 



d e (smddg{-)). 



(4) 



This equation is supplemented by the following con- 
ditions: 1) Conservation of volume, / sin^d^C = 0; 2) 
The position of the "contact" at 8 = 8 C is given by 
(R + ((8 c ))cosd c = R (see Fig. |2J}; 3) The contact an- 
gle is 180°, implying that ('(8 C ) = Rtnn9 c . 

These equations are readily solved and we get 
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1 + cos 8 
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1 — cos 8 C 



with 



2k 2 R 2 /3, 



(5) 



(6) 



as derived in [5J and checked experimentally in [5]). We 
then calculate the static sag of the top of the droplet 
e = — C(tt), and that of the center of mass eg, yielding 
respectively: 
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where we recover the logarithmic correction assumed in 
Eq. (JTJ). In order to interpret these equations, we invert 
Eq. [SI which yields to leading order: 




Fig. 2: Solide circle: Schematic profile of the sessile droplet 
compared to the shape of a sphere of same volume centered in O 
(dashed circle) . 8 C represents the location of the "contact" line 
of the drop with the surface. The sag of the drop is measured 
by the shift e of the top of the droplet compared to that of the 
unsquashed sphere. 
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g = fc eff e G , 



with 



47T<7 
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(9) 



(10) 



Eq. 1101 can be interpreted balance between the 
liquid weight of the droplet and a non-linear elastic force 
fccff^Gi with an effective stiffness k c ff. Remarkably, the 
stiffness vanishes for small deformations, meaning that 
the liquid spring becomes softer when the deformation is 
weaker. During small amplitude impacts or oscillations, 
we therefore expect the characteristic time of deformation 
to diverge as 



Building on this qualitative argument, we now give 
a more quantitative calculation of the dynamics of the 
droplet. Using the same argument as for Hertz contact 
[TJJ , the existence of a characteristic timescale longer than 
the oscillation period y pgR? ju of the Rayleigh oscilla- 
tions of a free droplet allows one to use the static response 
studied in the previous section to calculate its dynamical 
properties. In the weak deformation limit the pressure in 
the liquid is given in first approximation by 2a /R, and 
Newton's law can be written as 



Met 



Mg-^(t), 



(11) 



where M is the mass of the drop, and E(t) = irR 2 8 2 (t) 
is the "contact" area. The problem is closed by deducing 
from Eq. ((BJ) and ([RJ the relationship between eg and S 
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Fig. 3: Dynamics of a weakly deformed drop on a non wet- 
ting substrate. For all graphs, times are scaled in units of 
l/wo = \J pR 3 /a, velocities in units of ^/a/pR, and lengths 
in units of k _1 = \J a I pg. Top: Contact time during an im- 
pact. Black squares: drops of radius 0.4 mm; Red diamonds: 
drops of radius 0.6 mm. Dashed line : prediction of (|13|l with- 
out gravity; Black and red solid lines: prediction of Eq. (|13|) 
with gravity for drops of radii 0.4 mm and 0.6 mm. Bottom: 
Small amplitude oscillation period of a sessile drop: black di- 
amonds: resonance frequency of a mercury droplet; red trian- 
gles: free oscillation frequency of a coated water droplet; black 
line: asymptotic expansion (|15p for small radii. Note that for 
both graphs, theory is plotted without any adjustable param- 
eter and that its range of validity is restricted to the limit of 
small deformation V/Vo <C 1 and nR <C 1 respectively. 
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Using Eq. (TTTT) and (fT2]) . we first study the impact of 
the droplet bouncing off a non-wetting surface. Just before 
impact we have e = 0, and e = V, where V is the impact 
velocity. The contact duration T c is calculated readily by 
noting that the right-hand term in Eq. (|12p derives from 
a potential energy U(e) defined by the implicit relation: 



dU = M{g-2a^{e G )/R )de G . 

Integrating the associated conservation of energy then 
leads to 



2de 



o yJV 2 - 2U(e)/M 



(13) 



where e m is the maximum deformation of the drop, defined 
by U(e m ) = MV 2 /2. In Fig. [3]we compare the prediction 
of equation (fT3"]) with data of ref. [7]) obtained for droplets 
of radii R = 0.4 mm and R = 0.6 mm bouncing over 
a superhydrophobic surface. The fairly good agreement 
between experimental points and theory provides a first 
validation of our approximations. In particular the model 
explains why the contact time increases in the regime of 
small deformations (V small) where we may naively have 
expected a constant contact time. 

Moreover, Eq. (fTTj) and (fL2|) allow one to study the 
small amplitude oscillations of the drop around equilib- 
rium. In this pursuit, we linearize Eq. which can be 
written as eq = —ui 2 e g , with 



3(T 



2irpR 4 



de G 



(14) 



where the subscript indicates that the derivative is evalu- 
ated at equilibrium. Differentiating Eq. (fj"2"j). one obtains 
after some straightforward algebra, 



3cr 
~P~R* 
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(15) 



To test experimentally this prediction, we measured the 
free oscillation period of Hg droplets as well as water 
droplets on superhydrophobic surfaces. The results are 
displayed in[3]b) as a function of the drop radius. It is ob- 
served that the dimensionless period is not constant, but 
increases (by a factor 3) as the droplet radius decreases. 
Again the system behaves in a softer way when it is more 
weakly deformed. In this limit, Eq. [T5]comparcs favorably 
with the data, without any adjustable parameters, except 
in the limit kR > 1 for which the drop is strongly flat- 
tened by gravity and cannot be described within a weak 
deformation approximation. 

The later result contradicts an earlier work on the oscil- 
lations of sessile drops where no logarithmic behavior was 
found [14] . We argue that this discrepancy is due to an in- 
consistency in the calculation presented in [14] where the 
eigenmodes of the droplet were decomposed on a single 
spherical harmonic even though the system is no longer 
rotationally invariant in the presence of the substrate. We 
show here that reintroducing the full set of spherical har- 
monics allows us to recover Eq. (fTB"]) and extends it to 
higher order modes. 

For high frequency modes, the quasi-static approxima- 
tion performed above is no longer valid, and one needs 
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in principle to solve the full Rayleigh-like problem in the 
presence of boundary conditions imposed by the substrate. 
In |14j , it was proposed to replace the boundary condition 
by a pressure field P c (r2) located on the contact region. 
In the limit of small droplets, the contact area is small 
and the pressure field can be approximated by a Dirac 
distribution 



F F 

p c (n) = ^|<j(n) = ^^2Y i °(e)Y e °(e 



0)*, (16) 



where F c is the force exerted by the substrate on the 
droplet and where we have decomposed the angular Dirac 
5 on the basis of spherical harmonic Y/™ (by symmetry 
around the z axis, the to 7^ do not contribute). Us- 
ing the definition of the spherical harmonics, we have 
Y £ °(0 = 0) = v/(2£+ 1)/4tt, hence 



i=0 



-in 



(17) 



The presence of this pressure field modifies the Laplace 
condition, and we now have 



Pr 



(2 + A e>¥ ,) C = AP + pgR cos 0, 



(18) 



This equation mixes the effects of gravity and oscilla- 
tions. However we can decouple them using the linearity 
of Eq. 13 We thus set ( = ^°)+C (1) and F c = P c (0) +P c (1) , 
where the subscript (1) denotes the static (dynamic) 
part of the deformation. Let's now project this equation 
on the basis of the spherical harmonics. Since P c possesses 
non-zero components on all Y®, Q 1 ' must be written as 

C (1) (°) = Y,tCi~\t)Y?(n). This point is the main dif- 
ference with the analysis of [33] where a less general form 
C^ 1 ) = £i(t)Y?(n) was assumed, in contradiction with 
the decomposition of P c on the full set of I values. 

In the case of an inviscid fluid, the flow in the bulk can 
be derived from a velocity potential %p(r,t) solution of the 
Laplace equation Ai/j = 0. Using the boundary conditions 
at the surface of the droplet d r ip(r — R, f2) — dt((n), 
the decomposition of ijj of the spherical harmonic can be 
readily expressed in terms of the £(l)^'s and we have 



iR l - 



(19) 



Using the time-dependent Bernoulli equation pdtip = 
—p, we can express the pressure field in Eq. (|18[) in terms 
of Q. Projecting Eq. (| on the basis of the spherical har- 
monics and using the condition Ag^Y™ = —£(£ + 1)Y™, 
we get an expression for Q: 



•d) 



a £(£-!)(£ + 2) 



uj 2 /uj 2 



2£+l 
An 



(20) 



with 



a I pR? . To obtain the excitation spectrum of 



the droplet, we use the condition Q^'iQ = 0) = imposed 
by the substrate. With the expression of Q 1 obtained 
above, this yields the following equation for the eigenfre- 
quency u 



£(2£+l) 



^/wg + 2) 



0. 



(21) 



The sum in Eq. (|2ip is divergent and it must be reg- 
ularized. The origin of this divergence comes from the 
singular nature of the 8 pressure-field used to describe the 
contact between the substrate and the droplet. Indeed, 
modes with high values of £ oscillate very fast spatially 
(on angular ranges of the order 9 ~ 2n / £) and can there- 
fore be sensitive to the actual shape of the contact region. 
This suggests to introduce in the sum (|2ip a cutoff at 
^max — — {kR)^ 1 ! as done in the static case. Using 
this assumption, we reformulate Eq. (|21[) as 



E 



£{2£+l) 



lu 2 /uj%- £(£-!){£ + 2) 



= 0. 



(22) 



To get a further insight on the solutions of this equation, 
we recast it as 



E 



i) 



u 2 /ul~ £{£-!){£ + 2) £ 



1 



(23) 



The series in the lhs of the equation is now convergent, 
and we can safely take the limit £ max —> 00. As for the 
rhs, using the definition of the Euler constant 7, it can be 
expanded as Y^e=T = W^max) + 7 + o(l). Finally, we 
can formally rewrite Eq. (|2"3"]l as 



F{u 2 /ul) = 2 (ln(i n 
with F being defined by 



7), 



m - E 



£{2£+l) 



z-i(£ -!){£+ 2) £ 



(24) 



(25) 



Since £ max is large, the solutions of Eq. (1241) must be 
close to the poles of F located at the values zi = £(£—l)(£+ 
2) corresponding to the Rayleigh frequencies of a levitated 
droplet. Let us expand F(z) close to zi up to the constant 
terms. We have in this limit F(z) ~ £(2£ + l)/(z - z e ) + 
ci£ + o(l), where at is a number that can be calculated 
numerically {at = -1/6, a 2 = 2.0889, a 3 = 3.0927,...). To 
this approximation we get the following dispersion relation 
for the sessile droplet 



£(£-!)(£ + 2) + 



£(2£ + l) 
ln(^ ax e27-a,)- 



(26) 



Except for the case 1=1 for which z\ = 0, the eigen- 
frequencies of the sessile droplet are given by the Rayleigh 
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Fig. 4: Comparison between resonances of vibrated Mer- 
cury droplets from [14], and the Rayleigh prediction ue/uo — 
^/£(£ -!)(£ + 2). 



frequencies of a levitated droplet, up to a logarithmic cor- 
rection. The case £ = 1 corresponds to the low frequency 
mode studied in the first part of this Letter. Eq. ([To]) and 
(f^G)) coincide if £ max = v^e -1-7 / nR, which agrees with 
the general scaling introduced earlier to define £ ma x- 

We now compare this prediction to the measurements 
of ref. [2] . We argue that in this article the values of I 
were not attributed to the proper modes. Indeed, Courty 
et al. claimed that, by monitoring the center of mass of 
the droplet excited by the vibration of the substrate along 
the vertical direction, their detection scheme could only 
detect modes with even values of I and thus attributed 
the value £ = 2n to the n-th mode. On the contrary, we 
argue here that the apparatus could track all resonances, 
since all modes arc coupled to the £ = 1 spherical har- 
monic describing the motion of the center of mass. We 
thus attribute here the value t = n to the n-th resonance 
observed experimentally Using this new attribution, we 
observe a fairly good agreement between the position of 
the resonances observed by Courty et al. and the predic- 
tion of Eq. (gH) (see Fig. g]) 

The comparison between the predictions of Eq. (|lip and 
(112[) with existing experimental data support our model of 
logarithmic spring to describe the dynamics of a weakly 
deformed droplet. In future works, we plan to extend our 
findings to other experimental situations which are yet 
unresolved. Among these, we will include hysteresis of 
the contact angle observed on real surfaces, which might 
explain the critical velocity below which the drop does not 
bounce anymore [15J . We will also study the impact on 
clastic or liquid surfaces. This would for instance help 
understanding recent experiments on droplets walking at 
the surface of water (T5J. Finally, it could be interesting 
to extend the present analysis to higher deformations to 
probe impacts at high velocity and explain the saturation 
of the contact time observed in [3] . 
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